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Electronic current in a nanometer-size rod is theoretically investigated by an eigen-channel 
decomposition method in nonequilibrium Green's function formalism. Physical properties, such 
as the local density of electrons and local current, are decomposed into contributions of eigen- 
channels. We observe that the evanescent modes and nonlinear conductance are enhanced in 
d-orbital systems, and the structure of the transmission function, local current density, and 
penetration depth are discussed. The two effects of the ill-contact at electrodes in d-orbital 
systems, evanescent modes and the nonlinearity of conductance, are regarded as originating in 
the peak structure of the transmission function of eigen-channels in the energy region between 
chemical potentials of left and right lead wires. 
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1. Introduction 

Nonequilibrium current through atomic-scale contact 
is very important in nanoscicnce and nanotechnology, 
where we must know in detail the effects of electrodes. 
Electrodes themselves can be scatterers for electron 
transport even if atomic structures are perfectly ordered 
at electrodes. A decay mode is found in electronic current 
in nanometer-scale conductance, when the conductor has 
freedoms of scattering path at each layer. For example, 
a decay mode is investigated in a monosite chain model 
with multiple orbitals for each site.^ The decay mode, 
evanescent wave, is known also in optics as a nonpropa- 
gating mode without dissipation and is of great interest 
for a possible foundation of nanophotonics. It appears in 
a nanometer-scale region of material surface, when the 
light is totally reflected at the surface.^ 

Realistic nanometer-scale conductors, on the other 
hand, have been studied extensively for their unique 
structural and transport properties.'^ One can find, for 
example, quantized conductance in a monatomic wire 
through a ballistic channel and resonant conductance 
through a quantized state with strong scattering at elec- 
trodes. In nanometer-scale systems, an atomic struc- 
ture is crucially important for transport properties.'* The 
structural investigation is essentially important for dis- 
cussing quantitative properties of realistic conductance 
of nanoscale rods. For example, gold appears as a heli- 
cal multishell nanowire,^ whose formation was explained 
with a two-stage process in nanometer-scale phenom- 
ena.® 

We like to discuss, in this paper, the characteristic 
properties generic in several nanometer-scale rods. Par- 
ticularly, the possible effect of ill-contact at electrodes 
in nanometer-scale conductors is investigated, as a com- 
bined effect of structural and transport properties. The 
evanescent wave current appears, when we carry out 
a channel-decomposition analysis with nonequilibrium 
Green functions, and the nonlinear conductance under 



a finite bias voltage is observed. 

We will present, in this paper, calculations of an 
electronic structure, for d-orbital systems of rodhke 
nanowires of different lengths and widths, which gives 
prototypical cases for real materials. The generic ill- 
contact in the electrodes is enhanced in d-orbital systems 
and becomes an origin of evanescent waves and nonlin- 
ear conductance. This paper is organized as follows; In 
§ 2, we review briefly the formalism of the nonequilib- 
rium Green's function and present the eigen-channel de- 
composition of Green's function, local current and other 
physical quantities. We show several calculated results 
on the decaying properties of electron density and the 
nonlinearity oi I/V curves in § 3. Section 4 is devoted 
to the analysis and discussion of the origin of nonlin- 
ear conductance and evanescent modes. The penetration 
of evanescent modes is also discussed. The discussion is 
focused on two effects of ill-contact in electrodes, the 
evanescent modes and the nonlinear conductance. The 
conclusion is given in § 5. 

2. Theory of Conductance in Nonequilibrium 
Green's Function 

We proposed a method of the eigen-channel decom- 
position of several properties of quantum transport in 
nanometer-scale conductors by using the nonequilib- 
rium Green's function (NEGF) formalism.^'® The eigen- 
channel is the transmission eigenstate obtained by the 
unitary transformation of the original channels,^ whose 
concept may be useful for analyzing current through 
nanometer-scale metal conductors. 

Brandbyge and coworkers obtained the expression for 
the eigen-channel decomposition of the local density of 
states (LDOS).^*^' We generalize their formulation and 
provide the decomposition of NEGF and transport quan- 
tities, such as the density of electrons or local current, 
into contributions from respective eigen-channels. This 
formulation is also applied to the case of finite bias volt- 
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age. Information on the eigen-channel decomposition of 
focal current may be necessary for analyzing the vortex 
of current. Later in the present paper, we will clariiy 
the contribution of 'evanescent' modes by eigen-channel 
analysis. 

2.1 Nonequilibrium Green's function 

We use a nanometer-scale metallic rod whose entire 
system consists of three parts, the conductor (C), the 
left lead wire (L), and the right lead wire (R). The two 
lead wires, L and R, are semi-infinite and connect with 
the central conductor C. 

The retarded and advanced Green's functions and 
G"^ in the conductor C can be obtained within the NEGF 
theory as'^'^ 



is an eigen- vector of the matrix t]^LiRL- The incoming 
wavefunction from the right electrode, ur„, is defined in 
the same manner. 

We can define the decomposition of the lesser self- 
energy into contributions from respective eigen- 
channels ULn and ur„; 

S< = I/lFl + ifnVn 

n 

= E^Ln + E^R- (7) 

n n 

This definition is a basic idea and, then, the decomposi- 
tion of the lesser Green's function G< can be introduced 

as 



G-f-iE) = (E-Hc-i:^^^E)-i:riE)y , (1) a< = G-E<G- = ^G<„ + ^G<„, (8) 



(3) 



i {fL{E)Tj^{E) + fn{E)Tn{E)} , (4) ^^ere pL„(i, E) is defined as 



where Hq is the Hamiltonian matrix of the conductor C 
and its matrix dimension Nc is given as 

A^C = (the number of atoms in the conductor) 

x(the number of orbitals per atom). (2) 

The coupling of C with the L(R) lead wires is rep- 
resented by the retarded/advanced self-energies, E^'J^^, 
which have non-zero elements only near the electrodes. 

The lesser Green's function and the lesser self- 
energy E< are defined as ^'^ 

G<{E) = G^{E)^<{E)G^{E), 

S<(£) 

where ^h(R) is a semi-positive-definite matrix given as 
rL(R) = ~ ^L(R)^) and /l(r) is the Fermi-Dirac 

distribution function in the unperturbed lead wire L (R) 
of infinite size. The lesser self-energy represents the rate 
of electron scattering by the electrodes. 

We will calculate the sclf-cncrgics by an efficient re- 
cursion method^^ with an appropriate smearing factor 
i5 = 1.0 X 10~^Ry and the temperature of electrons in 
the electrodes is assumed zero. It should be noted that 
strictly localized states in C without interaction with 
incoming/outgoing modes in the lead wires are not in- 
cluded in G<}'^ 

2.2 Channel decomposition of physical properties into 
eigen- channels 

The channel-decomposition method was introduced for 
the local density of states (LDOS) by Brandbyge and 
coworkers.^''' Now, we arc generalizing the method by 
defining the channel decomposition of NEGF. The theory 
is formulated with real-space (atomic) bases. 

The transmission amplitude matrix from L to R, tjo,, 
and that from R to L, ^lr, are given as^^ 



where G^f^^-^^^ is defined as 



(9) 



By using this eigen-channel decomposition, the local 
electron density pi{Ey'^'^'^ is decomposed into eigen- 
channel components; 

Pi{E) = ^lmG<{i,i){E) 

PRn{i,E), (10) 



PLn{i,E) = — ImG<„(^,^)(i?) 



1 
1 

2^ 



/LRe[(G«r£«L„4„r£G^)ii] 



/LRe[(G'^r2uLn)^(G^r£uL„)j3ll) 



and pYin{i,E) is defined in the same manner. 

The local current ii^j(S)^'^'^^ is also decomposed, 

when the energy E lies within the chemical potentials 
at the left and right electrodes (/xr < E < pi^); 



.Go, 



ii^j{E) = 2^Re[Hcii,j)G<AE)] 



Go 



(12) 



tRL 
thR 



(5) 
(6) 



An incoming wavefunction of the cigcn-channcl "n" from 
the left electrode is denoted as the vector Uhn, which 



where 

in,i^j{E) = 2Re[HjiG<{i,j){E)] 

= -2Im[if,i(G^rJwLn)i(G^rjuL„)*]. 

(13) 

Here, Go is the quantum unit of conductance including 
spin degeneracy (Gq = ~ 0.775 x 10"'' for spin- 
unpolarized calculations.''*'^^ 

The eigen-channel decomposition of LDOS by Brand- 
byge et al.^^ is given by substituting /l = /r = 1 into 
the expression of the local electron density, eq. (10). 
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3. Evanescent Wave Modes cind Length Depen- 
dent Conductance 

The presence of the evanescent wave current and non- 
linear conductance is demonstrated, in this section, in a 
nanomctcr-scalc conductor with an electronic structure 
of s and/or d orbitals, which may be a prototypical model 
for realistic metals. Here, we adopt the Hamiltonian of 
spin-unpolarizcd bulk (fee) gold (Au) only with s-s or 
d-d hopping integrals without s-d hybridization in the 
framework of the tight-binding linear muffin-tin orbital 
(LMTO) method with the atomic sphere approximation 
(ASA),^^"^^ specifically the first-order Hamiltonian H^^\ 

3.1 Gold nanorods: calculated system,s 

Calculated systems are nanometer-scale conductors: 
(1) systems of semi-infinite size extending to ^; = -|-oo 
connecting with a left lead wire, as depicted in Fig. 1(a) 
and (2) those of finite size connecting with left and right 
lead wires, as shown in the inset in Fig. 2(b). The con- 
ductors and lead wires are quarried from an ideal gold fee 
lattice with a lattice constant of a = 7.755 a.u. Actual 
gold nanowires have helical multishell nanostructures^ 
whose formation process is explained with the two-stage 
process based on the quantum molecular dynamical sim- 
ulation.^ However, in the present paper, we use the above 
toy model of gold nanorods with an ideal fee structure, 
without lattice relaxation nor reconstruction. 

The left lead wire L consists of two layers of 64 atoms 
and the right lead wire R is the conductor itself of semi- 
infinite size (Fig. 1(a)) or the left and right lead wires 
L and R are the same two-layer blocks of 64 atoms of 
finite size (Fig. 2(b)). The cross sections of the wires C 
are 1x1,2x2, and 3 x 3 fee conventional unit cells. We 
denote the structures by C-1 x 1, C-2 x 2, and C-3 x 3, 
respectively. The lengths of the wires are 12 layers for 
C-1 X 1 and C-2 x 2, and 9 layers for C-3 x 3. Each layer 
consists of two sublayers because of the fee structure, one 
atom layer contains five atoms and the other one contains 
four atoms in each unit (1 x 1). 

The periodic boundary condition is imposed on the 
X ^ y plane with the primitive vectors (4a, 0, 0) and (0, 
4a, 0), and only the F point in the kx-ky plane is used 
for the electronic structure calculations. In other words, 
the cross section of the lead wires is of periodic 4x4. 

Effects of a bias voltage V are added to the Hamilto- 
nian. In the present work, the LDA potential is not de- 
termined sclf-consistently in conductors at a finite bias 
voltage and a uniform voltage drop of a bias potential is 
assumed within a conductor. The voltage drop in a con- 
ductor may not be uniform in actual nanometer-scale 
materials but may occur near the regions of contacts 
where charge distribution would change largely, and 
the self-consistent calculation for LDA potential in thick 
nanometer-scale wires may be important for conduct- 
ing detailed quantitative analysis. The structure of volt- 
age drop presumably depends on the conditions of elec- 
trodes and the conductor; for example, in ill-connected 
electrodes, the voltage drop occurs in the region of elec- 
trodes. 
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Fig. 1. (a) Calculated system of a C-1 x 1 structure. The con- 
ductor region (C) has k = 1,2,3....12 periodic layers. Here, a 

periodic layer is defined as two atomic layers (sec text for de- 
tails) (b) Decaying behavior of the averaged current mode, due 
to the presence of evanescent wave components. The number of 
electrons per periodic layer p'*^', defined in eq. (14), is plotted as 
a function of the number of periodic layers. No bias voltage is ap- 
plied, and electrons enter the system only from the left electrode 
(/l = 1, /r = 0). 



3.2 Decaying electron density in rods and evanescent 

mode 

We show the electron density for s- and d-channels in 
Fig. 1(b) for several semi- infinite size rods (Fig. 1(a)) of 
different cross sections, in which the local electron den- 
sity {pj (E)) is integrated with respect to energy and plot- 
ted as the sum of each periodic layer; 



/ + 00 
-oo 



jGfe— th layer 



(14) 



and we observe the decaying characteristics of electron 
density. Here, no bias voltage is applied, and electrons 
enter the system only from the left electrode (/l = 1, 
/r = 0). The quantity of p^'^^ shows the averaged char- 
acteristic among all existing modes. We also calculated 
the electron density profile in similar periodic systems 
with a larger cross section of the lead wires of 6 x 6, and 
no appreciable change was observed. Therefore, we can 
conclude that the above dependence on the cross sec- 
tion of conductors is universal but not a size effect of 
the adopted systems. The decaying characteristics of the 
electron density can be expressed as 



P 



(k) 



P 



[oo) 



p(o)e- 



(15) 



and we can define the penetration length ^. 

Figure 1(b) indicates the existence of components of 
evanescent wavefunction among the modes, which are 
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Fig. 2. Length dependence of transport properties for the d- 
channel system of C-1 x 1 cross section: (a) Transmissions T{E) 
of the six-layer conductor (the broken hne) and the 1-layer con- 
ductor (the solid Hnc) at = V, both assuming /ir = —6.6 
cV. The transmission of tlic sliortcr rod is larger, owing to tiie 
contribution by evanescent wave modes, (b) I-V curves of the 
six-layer conductor (the broken line) and the two-layer conduc- 
tor (the solid line). Here, we assume = —6.6 cV. The longer 
rod shows a higher nonlinearity because of the enhanced spiky 
structure of the transmission function. The inset shows the ge- 
ometry of the six-layer conductor. 



ncar-ficld standing wavefunctions decaying exponentially 
with increasing distance from the electrode. The compo- 
nents of evanescent wave modes decay owing to backscat- 
tering and. therefore, only Bloch wave modes reach the 
right end {z — -|-oo). In other words, when the conductor 
length is comparable to or smaller than the decay length, 
the evanescent wavefunctions can reach the opposite elec- 
trode. Therefore, one can expect the length dependence 
on conductance. Note that the comparison between the 
s- and d-orbital channels is given later in § 4.3. 

3.3 Nonlinear conductance with finite bias voltage 
The transmission function T{E) is defined as^'^'^^ 

T{E) = ii[G^{E)T^{E)G\E)T^{E)\. (16) 

By using eqs. (5) and (6), the transmission function T{E) 
can be decomposed into the eigen-transmission |t„P with 
the unitary matrix C/l in terms of diagonalization of 



T{E) 



tr 



El 



= diag{|T„|^}. 



(17) 



The incoming wavefunctions of the eigen-channel from L 
(the left electrode) can be given by the components of 
the unitary matrix J7l or mlh- 

The current can be evaluated by using the transmis- 
sion function with the relation^' ^'-'^^ 

nV) = — / T(E)dE, (18) 

where V is the bias voltage applied to the system, /xl(r) 
is the chemical potential of the L (R) lead wire (/^l = 
liB. + eV,V> 0). 

The length dependence in the conductance of a d- 

orbital channel is shown in Figs. 2(a) and 2(b) for sys- 
tems of different rod lengths with a C-1 x 1 cross section. 
The calculated conductance in Fig. 2(a) shows that the 
shorter rod has a larger conductance owing to the pres- 
ence of the evanescent wave modes. In longer systems, 
the Bloch modes only contribute to conductance. Here, 
we notice a larger nonlinearity in model systems in the 
I-V curve of Fig. 2(b), because a strong resonance scat- 
tering in the d-orbital channel causes a strong energy de- 
pendence in the transmission function T{E), which will 
be seen in the next section. 

4. Origin of Nonlinecir Conductance and Evanes- 
cent Modes 

4-1 Quantization of one- dimensional modes along rod 
length in d-orbital channels 

Now, we investigate characteristics of the eigen- 
channels by generalized channel analysis in § 2. 

Figure 3(a) shows the cigcn-transmissions |r„p of the 
six-layer conductor in Fig. 2, which shows a strong energy 
dependence or a spiky structure of an energy interval of 
0.01 ~ 0.02 eV. The eigen-transmission within a narrow 
energy range (—5.8 eV < E < —5.4 eV) is depicted in 
Fig. 3(b), in which we focus on the eigen-channel 'A'. The 
oscillatory behavior of the transmission ftmction of the 
eigen-channel A is due to the resonance and off-resonance 
of one-dimensional Bloch modes. 

The d-bands in bulk systems have a width of about 
4 eV. Then, the finiteness of the rod radius in the x- 
y plane makes each band split. This splitting causes a 
complicated band-crossing since some d-bands exhibit 
the dispersion relation dE/dkz > 0, whereas others ex- 
hibit the dispersion relation dE/dk^ < 0. Then, a com- 
plicated branching occurs in d-orbital systems. We plot, 
in Fig. 3(c), the E-kz curve of the 1x1 infinite conduc- 
tor representing the branch of a one-dimensional mode 
'A' corresponding to the 'A'-channel. 

The right vortical axis refers to the phase A'cfffl x kz 
of the branch 'A' in an infinite periodic wire of fitted 
effective periodicity length A'cff' = 4.88. In fact, the peri- 
odic unit length A'effO of 4.88 is determined so that the 
perfect transmission (|t„P = 1) occurs at the phases of 
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Fig. 3. Eigen-transmission for the d-channel system in a narrow 
energy range, (a) Eigen-transmission of the six-layer conductor 
at y = V and (b) eigen-transmission as a function of energy 
near the d-band bottom. The eigen-channel A is shown by a thick 
line, (c) E-k curve of an infinite conductor and the phase of one- 
dimensional modes (right-hand side) in the conductor, which is 
calculated as N^gaxkz (N^g = 4.88), where N^ga is the effective 
length of the conductor. The open circles in (c) are the points 
used in the fitting procedure of A'off and correspond to A2, A6, 
and A7 in (b). 



S — TT, 2tt, and Stt. This N^ff value (4.88) is nearly iden- 
tical to that of the conductor part {N = 6). The fitted 
length unit of A^off — 4.88 is actually shorter than that of 
the actual conductor part {N ~ 6), since the wavefunc- 
tions are scattered in the regions near the electrode and 
the resultant effective conductor length should be short- 
ened from that of the ideal one. The above analysis in- 
dicates that each eigen-channel (branch) of the d-orbital 
system has peaks of the transmission corresponding to 
a phase at iVoffO x = mr. The energy interval be- 
tween branches is determined by the size of the rod and 
the structure of d-bands and should be inversely propor- 
tional to the radius of the rod. In other words, we have 
observed two types of nanometer-scale rod quantization; 




L2 L1 1 2 3 4 5 6 R1 R2 
Index of layer 

Fig. 4. Eigen-channel decomposed density of electrons for the 
eigen-channels from Al to A7 in Fig. 3(b). Electrons enter the 
system only from the left lead wire (/l = 1, /r = 0). Broken 
lines show the function a sin^ ( -j^^ (x — 3.5 + )) -I- b with the 
layer index x and n = 1,2, and 3 for A2, A6, and A7. The 
adjustable parameters (a,b) are (8.0, 0.9), (3.9, 1.1), and (5.6, 
1.2) for A2, A6, and A7, respectively. These lines suggest that 
the charge densities at A2, A6, and A7 may have one, two, and 
three peaks, respectively, as the phase in Fig. 3(c) shows. 



the quantization of the radius of the rod (branch) and 
that of the length of the rod (peak structure of eigen- 
transmission) . 

The oscillatory behavior of resonance and off- 
resonance in s-orbital channels is relatively weak;^^ one 
branch extents over a much larger energy range. For ex- 
ample, the energy width of one branch is 2 '--^ 4 eV and 
the energy width of the s-band is about 10 eV in Au. The 
energy interval of transmission peaks is wider, and the 
difference between the maximum and minimum values is 
small. This oscillatory behavior difference may be con- 
sistent with the decaying profile of electron density in s- 
and d-orbital channels in Fig. 1(a), i.e., the evanescent 
wavefunctions appear in order to connect with wavefunc- 
tions in the conductor and lead wires in the electrode 
region. In the d-orbital channels, more evanescent wave- 
functions are necessary than in s-orbital channels, be- 
cause the wavefunctions have more freedoms of orbital. 
The strong resonance scatterings causing the large oscil- 
latory behavior in d-orbital channels are consistent with 
previous works on transition metals. ^^"^"^ 

The evanescent wavefunctions of d-orbital channels 
penetrate much deeper and are more scattered at the 
electrodes, in order to connect with Bloch modes. This 
we call the 'ill-contact' of d-orbital channels. The channel 
decomposition of the local density of electrons is shown 
in Fig. 4 for several energy points Al~A7 of the eigen- 
channel 'A' in Fig. 3(b). We can see the characteristics 
consistent with the above mode analysis results; e.g., the 
point Al is placed at the tail position of the branch 'A' 
in Fig. 3 (b) and its local density of electrons exhibits a 
typical decaying (evanescent) behavior. The point A2 is 
placed at a peak of the eigen-transmission of the branch 
'A' in Fig. 3(b) that exhibits a nodeless standing wave, 
which is consistent with the assignment of the phase of 
S — TT in Fig. 3(c). The density profile shows that the 
connection of transmission in C and R or L is poor and 
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Fig. 5. Voltage dependence of transmission function T(E) for 
three cases: Vi = 0.291 V, V2 = 0.332 V, and V3 = 0.367 V. 
The chemical potential fi^i = —6.6 eV and respective energies 
/JR + eV are shown by broken lines. The transmission function 

T(E) between and = MR + contributes to the current 
for the bias voltage V . The peak shown by a broken circle creates 
the structure of ^j^^^ in the I-V curve. The peak structures 

marked by solid circles contribute strongly to (jf)^^,' 



does not smear out from C. 

4-2 Nonlinear conductance and bias voltage dependence 
of transmission 
Let us discuss the finite voltage case. The transmis- 
sion function T{E) strongly depends on the bias volt- 
age V, as shown in Fig. 5. Furthermore, there exist two 
characteristic features of d-orbital channels: (i) the neg- 
ative differential conductance (3^ < 0) and (ii) the peak 
structure in the differential conductance. The differential 
conductance is separated into two parts: 

dI{V) 



GoT{E = MR + eV, V) 
Go f^'''+^^ dT{E,V) 



+ 



/ 



dJ 
dV 



+ 



dl_ 



dV 



AT 



dE 



(19) 



The first term (^)j, is the contribution from a new 
structure of the transmission function T{E) entering the 
energy window (/xr + eV > E > /zr) from the higher 
energy side and is always positive. The second term 
(57) AT ^® contribution from the change of the over- 
all structure of the transmission function T{E) in the 
energy window and can be negative. Therefore, the neg- 
ative differential conductance is due to the strong voltage 
dependence of the overall structure of the transmission 
function, which is caused by the channels confining the 
local electron density within the conductor. 

The sensitivity of the resonance to the bias voltage is 
due to the small width of branches of d-bands as a result 




10 50 

Number of atoms on the cross section A'atom 

Fig. 6. Dependence of penetration lengths ^ on the number of 
atoms on cross section of the conductor iVatom for d-orbital chan- 
nels for the local electron density p^^^ in Fig. 1(b). The points 
in the figure indicate the values calculated using eq. (24) with 
error bars and the dotted line shows the behavior of -W^tom^" 
determined by the least-square fitting procedure. 



of ill-contact of d-orbital channels. The peak entering the 
energy window of ij,r + eV > E > fiR from the higher 
energy side gives rise to the peak of Therefore, the 
peak structure of differential conductance originates from 
that of the transmission due to the resonance scattering. 

4-3 Quantization d,ue to finiteness of radius of rod 

The mode analysis in § 4.1 indicates that the decaying 
behavior in Fig. 1(b) appears not only in specific decay 
modes but also in the sum of various modes, such as 
those in Fig. 4. 

Now, we explain the physical picture of the difference 
between s- and d-channels for the decaying behavior in 
Fig. 1. Here, we denote the number of bases of elec- 
tron wavefunctions on the cross section of the conduc- 
tor iVcross, the energy width of the s- or d-band T4^band 
and that of one branch Wbranch- For example, one can 
see from Fig. 3 that Wband — 4 eV (the width of the 
d-bands) and T4^bianch — 0.2 cV (the branch A). iVcross 
should be proportional to the number of atoms on the 
cross section of the conductor A^atom- 

A mode can contribute to the current, when it connects 
with a Bloch wavefunction in the conductor C in the 
contact region. The number of one-dimensional bands in 
the conductor C can be estimated as 

^^branch 



If the ratio 



R: 



Wband 



Wbranch 



(20) 



(21) 



W^band 

is small, it is difficult for incoming modes to connect with 
the one-dimensional Bloch modes in C and majority of 
the incoming waves may be backscattered. The ratio R 
is, in our examples, of the order of 0.2/4 ~ 0.05 in the 
d-orbital channels. 

The above explanation is consistent with the data in 
Fig. 1, since the contribution of evanescent modes can be 
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estimated as the difference in local electron density be- 
tween the left end (z = 1) and the right end (z = oo). In 
the case of the C-3 x 3-d system, for example, the num- 
bers of electrons at the left and right ends are nioft ~ 350 
and bright ~ 180, respectively. The contributions of (one- 
dimensional) Bloch modes (nBioch ~ "-right) and evanes- 
cent modes (nEvan ~ f^ioft — "right) are almost the same. 
The evanescent modes are very important for the d- 
orbital channels, since nEvan/^-Bioch ~ 1- 

The same analysis of the s-band system indicates that 
the evanescent modes are less important ("Evan/"Bioch < 
1) than that in the d-orbital channels, which is consistent 
with the fact that the ratio R — M^bianch/W^band is larger 
than that in the d-orbital channels as i? ~ 0.2 ~ 0.4 
because Wband — 10 eV and IVbranch ~ 2 ^ 4 eV. 

The discussion of the connection is supported by the 
dependence of penetration length on the number of 
atoms in the cross section. For simplicity, an extreme 
case with dispersionless sub-bands is considered for d- 
orbital channels. In this case, the mean energy interval 
-Eint at Fermi energy is scaled by 

Eint cx — (22) 

across 

Because the penetration length ^ can be estimated as 

1 



int 



we obtain 



oc 



(23) 



(24) 



Figure 6 shows the dependence of the penetration length 
^ on the number of atoms in the cross section A^atom for 
d-orbital channels for the local electron density p^'^^ in 
Fig. 1(b) with eq. (15). Because the penetration length 
is almost scaled by eq. (24), we can conclude that our 
simple explanation of the connection captures an essen- 
tial nature of evanescent waves. 

It should be noted that a discussion of scattered waves 
gives the same dependence of the penetration length. 
Given incident waves are backscattered at edges of the 
conductor, the collision rate per periodic layer is propor- 
tional to 1/L, where L is the radius of the conductor. 
Because the penetration length is proportional to the in- 
verse of the collision rate and Naoss cx L^, eq. (24) is 
obtained again. This supports our discussion in this sub- 
section. 

4-4 Eigen- channel decomposition of local current and 
evanescent wave modes 

The backward flow cannot be observed even if we 
would try to observe the total current flow at a certain 
layer position because of the conservation law of current. 
In other words, the net flow conserves at each atomic 
position. Furthermore, no backward flow is observed in 
the total current flow because of strong forward com- 
ponents of the other modes. Therefore, it is essential to 
decompose the local current into the eigen-channel com- 
ponents and, then, the eigen-channel current can show 
the backflow related to the evanescent modes. 

We show that the local current defined by eq. (12), us- 
ing the generalized channel analysis, gives the backward 




c 



L2 



D 



Fig. 7. Channel-decomposed local current at the energy point 
Al near the electrode region between the conductor and the 
load wire; (a) Geometrical structure near the electrode, (b) The 
channel-decomposed local currents are plotted as arrows and the 
numbers indicate their values. The lengths of arrows indicate 
their amplitudes on the logarithmic scale. The backward current 
flows from atom D to atom E (2 X 10~^) and from atom B to 
atom D (4 X 10-^) are several orders of magnitude weaker than 
the forward flows. 



current flow. Figure 7 shows the channel-decomposed lo- 
cal current in eq. (13) for the channel Al in Fig. 3. The 
backward current flow is found, for example, between 
the atoms B and D and between the atoms D and E. 
The same analysis for the channels A2 and A3 (figure 
not shown) indicates that the backward current flow ap- 
pears commonly between the atoms D and E among the 
above three channels. The above result supports the dis- 
cussion in the previous section, in which the evanescent 
wave component should be attributed not only to de- 
cay channels, such as the channel Al, but also to many 
other channels. Moreover, the observed backward flow 
should be a net flow of the forward and backward cur- 
rent components and the backward components should 
be inherent among other regions. 

We should emphasize that such backward flow can be 
obtained, only with the channel decomposition method 
given in the present paper. 

5. Conclusion 

An eigen-channel analysis method has been extended 
to general transport quantities based on the nonequi- 
librium Green's function formalism and we studied the 
transport properties of rodlike nanometer-scale metal 
wires from the viewpoint of "ill-contact" in d-orbital 
channels. The ill-contact generates two different charac- 
teristics in transport phenomena of d-orbital channels; 
the evanescent wave modes and nonlinear conductance 
at a finite bias voltage. 
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Although the present transport calculation has been 
carried out only in cases with model structures, the chan- 
nel decomposition method is general and applicable to re- 
alistic materials whose atomic structure is determined by 
an energy functional with an electronic structure. Such 
combined simulations of structural and transport proper- 
ties during the formation of nanoscalc struc;tiircs should 
be a foundation of nanoelectronics, together with com- 
bined experiments of structural and transport observa- 
tions. 
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